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ABSTRACT 

APPLICATION OF THE FAST FOURIER TRANSFORM TO THE 
COMPUTATION OF ELECTRIC FIELDS AND ASSOCIATED 
HEATING PATTERNS IN HUMAN TISSUE 

by 

i 

EDWARD MICHEAL SCHAEFER 

Chairman of Supervisory Committee: David L. Johnson 

Department of Electrical 
Engineering 

ReLative heating in human tissue is defined in terms of the eLectric 
fields. Analytical expressions for the electric fields in planar layers of 
human tissue due to a known aperture distribution at the surface of the 
tissue are derived. The Cooley-Tukey algorithm for efficient computation 
of discrete Fourier transforms is illustrated for a one-dimensional 
problem. Flow charts for the use of the fast Fourier transform (FFT) 
in relative heating problems are developed. Results from computation 
of relative heating by the FFT are compared to similar results obtained 
by the Guassian quadrature transformation procedure. 
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CHAPTER 1 



INTRODUCTION 

1, 1 Our Changing Environment 

More technological advances have been made in the Twentieth 
Century than in all of man's past history. These advances have affected 
nearly every aspect of existence on this planet as weLl as fulfilling the 
age-old dream of setting foot on other heavenly bodies, Man, as a result, 
enjoys a longer and easier life. 

These phenomenal advances have not been made without some 
potentially harmful side effects to the environment. Factories have con- 
tributed to air and water pollution. Questions are being raised about the 
relative merits of supersonic travel versus the adverse effects of the 
resultant sonic boom. Every man has a car in his garage but the country- 
side is being ruined by ribbons of concrete. Many recent discussions have 
centered on the adverse effects of electromagnetic and particulate radiation 
on biological systems in general and human tissue in particular, but there 
are also many advantageous uses for electromagnetic and particulate 
radiation such as the use of electromagnetics in diathermy treatment of 
many musculoskeletal disorders and the use of nuclear energy to generate 
electric energy. 

1. 2 Heat for Therapy 

Man has searched for more effective methods of applying heat to 
his body for various aches and pains since he first discovered what the 



2 



heat of the sun could do for him. This search has led man to investigate 

nearly every method of heat production including the sun, hot water, 

infrared radiation, ultrasonics and electromagnetic radiation. In recent 

2 

years, there have been many reports of successful applications of electro- 
magnetic heating of deep body tissue for therapeutic purposes (diathermy). 

Since diathermy has proven to be effective for treatment of many 
musculoskeltetal disorders, the search is continuing in the Department 

3 4 c 

of Physical Medicine and Rehabilitation at the University of Washington * 9 

in an effort to find the most efficient method of realizing deep tissue heating 
without damaging the skin and subcutaneous fat. 

Computer solutions of the mathematical models of the diathermy 
problem have in the past been confined to investigating only a rather 
limited set of points in the volume of the models. This research will 
investigate an algorithm which can more efficiently provide a much more 
complete picture of the heating. 

A direct consequence of diathermy research will be the develop- 
ment of the analytic methods which can be used in the study of the adverse 
effects of electromagnetic radiation and as a quantitative aid in the estab- 
lishment of radiation safety standards. 

1. 3 Radiation Safety Standards 

The relative importance of determining the biological effects of 
electromagnetic radiation can be ascertained by the attention given to 
the problem in recent years: 
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a. Numerous congressional hearings were held on the subject in 
1967 and 1968. 6> 7 

Q 

b. The Radiation Control for Health and Safety Act of 1968 
(RCHSA-1968) was enacted as a result of the congressional 
hearings . 

c. The Armed Forces have, since 1957, held an annual Tri- 

Service Conference on the Biological effects of Microwave 

9 

Radiating Equipments. 

d. There is continuing research into the possible uses of electro- 
magnetic radiation in the field of medicine. 

e. There has been a great deal of research conducted on the use 

of electromagnetic heating of deep body tissue for therapeutic 

, ,. x 3, 4, 5, 10 

purposes (diathermy). 

Discussion and research into the biological effects of electromagnetic 
radiation have resulted in both the United States and the USSR accepting 
radiation safety standards for prolonged exposure to radiation. 

1. 4 The Need for Research 

That there is a need for further research can be evidenced by the 
fact that the safety standards set by the two countries differ by three orders 

g 

of magnitude. In fact, the RCHSA-1968 directs the Secretary of the 
Department of Health, Education and Welfare to conduct studies of the 
problem and report the findings to Congress. 

How should the research be conducted? In the past, most of the 
research on the biological effects of electromagnetic radiation has been 
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qualitative rather than quantitative. Considerable research has been 
conducted on animals and the results have been extrapolated to man with 
a fair amount of success. But in fact, the very question of qualitative 
versus quantitative research has been debated at great lengths while the 
proponents of each side continued to conduct the research as best they 
could in view of the facilities at their disposal. 

1, 5 An Engineering Approach 

A reasonable engineering approach to the problem would incorporate 
both quantitative and qualitative research in a team effort by engineers and 
physiologists where members of each discipline would investigate that 
aspect of the problem to which they are best suited. 

The engineers are well suited to investigate the quantitative aspects 
of the problem by first constructing planar and cylindrical layered models 
which exhibit the electromagnetic properties of human tissue. This analy- 
sis would mainly encompass the problem of determining by analytical 
methods the electric and magnetic fields within the tissue layers of the 
models for different internal or external source configurations. 

When the models have been developed to the point where they reliably 
agree with theoretical calculations, the same modelling materials can be 
used to model the complex geometries of the human body. The answers 
derived from the models can then be correlated with data from qualitative 



methods. 
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1. 6 Research Objectives 

The approach outlined in the preceding section is currently in use 
at the University of Washington. The problem under investigation is the 
determination of the optimum frequency of electromagnetic radiation and 
applicator configuration to be used for therapeutic heating of deep tissues. 

The subject of this particular research is a comparison of the 
Guassian quadrature and fast Fourier transform (FFT) algorithms used 
to compute the electric fields induced in planar layers of human tissue. 
The vehicle problem is the computation of the relative heating pattern 
arising from a direct contact diathermy applicator. The applicator has 
a TE waveguide mode radiating into a bi-layered, semi- infinite lossy 
medium (See Figure 1. 1). A subcutaneous fat layer overlying a semi- 
infinite muscle is modelled by this medium. 

Current analysis utilizes the Gaussian quadrature method to per- 
form the necessary transformations. The economics of computation 
time in the continuing research program dictate that only a relatively 
few points in any plane of the tissue can be computed and those must be 
along an axis of symmetry. The fast Fourier transform will compute 
an entire array of points in any desired plane, thereby presenting a more 
complete picture, at a considerable savings on a cost per point basis. 
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Figure 1 . 1 

Model of Human Tissue 



1. 7 Content of Succeeding Chapters 

In Chapter 2 the electric field equations and relative heating 
expressions to be solved are developed from the literature. Chapter 3 
discusses the fast Fourier transform which will be used to accomplish 
the inverse transforms required to compute the electric fields for the 
relative heating patterns. 

The use of the FFT in this and another related problem will be 
presented by discussion and flow diagrams in Chapter 4. 

This research project analyzes the merits of using the FFT to 
compute the electric fields necessary to obtain the relative heating patterns. 
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The results from the fast Fourier transformations at several frequencies 
are compared in Chapter 5 to results obtained from Gaussian quadrature 
transformations already accomplished by Guy and McDougall. 

A brief discussion of the conclusions to be drawn from this 
research is presented in Chapter 6. 



CHA PTER 2 



RELATIVE HEATING IN PLANAR LAYERS OF HUMAN TISSUE 
2. 1 Objectives 

Relative heating from a rectangular aperture in an infinite conducting 
plane radiating into a two layered model of human tissue is defined in Sec- 
tion 2. 2 (See Figure 2. 1). The human tissues modelled are a subcutaneous 
fat layer of depth z ^ overlying a semi- infinite muscle layer. The succeed- 
ing sections develop the necessary electric field equations required for 
computation of relative heating. 

The development of the field expressions follows the approach of 
12 13 

Villeneuve and Guy. Villeneuve utilizes the Fourier transform 
approach in computing the effective input impedance of an open-ended 
waveguide radiating into a plasma environment. This problem arises in 
designing communication antennas on high speed reentry vehicles used in 
space exploration. 

Villeneuve, in his work, develops general expressions for the fields 
exterior to the waveguide; then, he looks at the expressions only at the 
aperture. Guy extends the expressions of Villeneuve by matching boundary 
conditions and derives the field expressions necessary to compute relative 
heating . 

2. 2 Relative Heating 

The criteria for the design of an optimum diathermy applicator 
for treatment of musculoskeletal disorders is that the heating per unit 
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Electromagnetic Source 




FIGURE 2. 1 



Tissue Model and Coordinate Definition 



10 



volume in the subcutaneous fat be significantly lower than the heating per 

unit volume in the musculature. The term relative heating arises from 

this criteria. Relative heating at any point (x, y, z) is defined as the 

heating per unit volume at the point (x, y, z) relative to the heating per 

unit volume at the point (0,0, z^), where z^ is the thickness of the sub- 

3 

cutaneous fat layer. Guy has shown that the heating per unit volume at 
the point (x, y, z) is given by 



= W + W + W 2. 2. 1 

T x y z 




W =?|Ep 2. 2. 3 

y 2 y 

w 2 ■ f P.f 2 - 2 - 4 



where 0" is the conductivity of the medium and the E*s are the electric 
field strengths at the point (x, y,z). From these expressions the relative 
heating at a point (x, y, z) in the fat layer is given by 



W (x, y, z) 

f 

W T (0,0, Zj) 

m 



2. 2. 5 



$ '(IeJ + fcf + pf) 






2 . 2 . 6 



and in the muscle by 



W T (x, y, z) 



m 



T W (0,0, z ) 

m T 1 

m 



2. 2. 7 



1 1 



p 




2 . 2 . 8 



where 



E = E (0,0, z ) 
xO x 1 



2 . 2 . 9 



E = E (0,0, z ) 
yO y 1 



2 . 2 . 10 




2 . 2 . 11 



The W (0, 0, z ) which appears in the p 
T i 



and p m 

f T 



express ions 



m 



normalizes the relative heating to unity at the fat-muscle interface. 

2. 3 General Solution of Electric Fields in Layered Media 

In rectangular coordinates, general solutions of the scalar wave 
equation may be obtained as a linear superposition of plane wave solutions.* 
This process leads to the development of a complete set of Fourier trans- 
form pairs which completely specify ail of the field components throughout 
the volume of the fat and muscle. 

The electromagnetic fields may be derived from two parallel 
vector functions whose direction in space is fixed. Assume these func- 
tions to be 



A = U i|; 



2. 3. 1 



x m 



F = U 



2. 3. 2 



x e 



where ib and Jj are solutions of the scalar wave equation 
e m 
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2 2 

[ V + k ] i|j = 0. 2.3.3 

A Linear superposition of plane wave solutions apropos to this 
problem (See Figure 2. 1) is, for Region (1) 



00 






1 



jW z -jW z .. 

J 1 J 1 j(ux + vv 

_ . . e + f, e )e J ' y; dudv 2. 3. 4 

m .2 J J 1 1 

(2tt ) -oo 



I 1 4 



00 



= r~2 



(2tt ) - 00 



f C 4 jW l Z - " jW l Z , j(ux + vy) 

J J (g j e + gj e )e dudv 



2. 3. 5 



And for Region (2) 



oo 



, 1 Hr, v j(ux + vy) ^2 Z J 

= 2 J J f 2 (u>v) e J e dudv 



(2tt ) -oo 



2. 3. 6 



oo 



jWz 



, 1 f f / \ j(ux + vy) J '2‘ , 

= s> J J g 0 (u, v) e JV e dudv 



2 J J & 2 
(2tt ) -oo 



2. 3. 7 



where 



/~ 1<2 2 ? 

W = /-s7 k - u - v 



2. 3. 8 



❖ 

k = 2 tt 



/ *1 
V l 1 € 



2. 3. 9 



From Maxwell's equations 



E = — — V xVx A - Vx F 
j0)€ 



2. 3. 10 



H - . Vx Vx F + Vx A. 
jcop. 



2. 3. 11 
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Making the appropriate substitutions and carrying out the indicated 
operations yields the following equations for the electric fields (See 
Figure 2.1). Region (1) 



oo 2 2 

1 r r W ! + v ^ jW n z -jW_ Z 

E. , U-4 (ft e 1 + f" e 1 >e> ( " x+v ’ r) dudv2. 

fx ,2 J jcoe 1 1 



(2tt) -oo JW 1 



oo 



E 



r c ~ uv f + - -j w i* 

jj— t (f i e tf i e * 



£y ( 2„) 2 1 

■* W 1 Z - ^1 Z . j(ux + vy) dudv 

-jW^gl g l n 



oo 



■uW, 



fz ,2 J J jwf, 1 

(2tt ) -oo J 1 



jW,z _ -jW,z 

1 






iW z -jW z 

,,+ J l - J 1 j(ux + vy) 

-J v (g j e + gj e ] e J dudv 



Region (2) 



oo 2 2 

* fe £W 2 Z ej (ux + vy, dudv 



mx 2 J J ico 2 

(2tt ) -oo J 2 



oo 



i W z 

2 j(ux + vy) 



my ( 2 , ) 2 -« j “2 2 
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oc - uW 



jW z .. 



Il ( — 2 ^ + jg ^ ,e 2 ej<UX+Vy>dUdV 2 



mz , 0 2 J J J ~2 

(2ir ) - 00 



3. 12 



3. 13 



3. 14 



3. 15 



3. 16 



3. 17 
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Matching boundary conditions at z = 0 and z = and solving 

the six resulting simultaneous equations gives the various f 1 s and g's. 

13 

This has been done by Guy with the following results for the fat-muscle 
case under investigation. 



Region (1), in the fat 
00 

E = \ f +,(u,v,z) e j|ux + "''’dudv 

fx kJ kJ l 

- OC 
00 
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OC 
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-oo 



2. 3. 18 



2.3.19 



2. 3. 20 



Region (2), in the muscle 



00 
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where 

4^ (u, v » z ) 



F (u, v, z) ^ 
F (u, v, 0) 



2. 3. 21 



2 . 3. 22 



2. 3. 23 



2. 3. 24 



-uv C 



[ F (u, v, z) - G (u, V, z)] 



1 5 



^(u, v, z) - 2 2 



(W, + v )F (u, v, 0) 



\ (U ’ V,Z)= ....2 2 



W 
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1 
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^(u,v,z) = T 5 e 

F (u, v, 0) 



-uve 



jW 2 (z - z L ) 



<l> 2 lu. v > z ) = z 2 



(W 2 + v )F (u, v, 0) 



[ F (u, v, z ^ ) + G (u, v, z ^ ) ] 



j W 2 (z - z L ) 
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X (u V z)=-^2 

2 V ’ ’ ' 2 2 



[■—■ F ’(u, v, z ) + G + (u, v,z )] 
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The F's and G ! s are given by 

jW,(z - z ) 

± + + J V i 

F (u, v, z) = (A sinW^z^ + j B cos W^) e 

-jW^z - Z] ) 

±(A sin W z^ + j B cos W^) e 



± + + + J W L Z 1 J W 1< Z - z l ) 

G (u,v,z) = [(A +B )sinW^z^ + jC e ] e 



-jW z -jW (z-z ) 
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+ 2 '!' ? 2 7222 

A = W t W 2 k « 2 T u v ± (WT + u )(W 2 + v ) 



2. 3.25 



2. 3. 26 



2. 3. 27 



2. 3. 28 



2. 3. 29 



2. 3. 30 



2. 3. 21 



2. 3. 32 
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B* = (W 2 + v 2 )(W^ + u 2 ) - u 2 v 2 ± W x W 2 k 2 t* 2. 3. 33 

± 2 * 22 * -> i i a 

C = WW k ejiWjk € x . 2.3.34 

Since the source field is at z = 0, then for alL z > 0 the fields 
must represent waves travelling and attenuating in the z-direction. That 
is, the branches of the propagation constants and W 2 must be 

chosen to yield these characteristics. Such a choice is 



W 1 


/ 2 . * 2 2 1 
= -/N/k - \1 - V 


0 < u 2 + v 2 < Re(k 2 ^') 
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where 
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2. 3. 40 
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2 IT f 




2. 3. 41 



c 



Now, the ? which appears in each of the electric field equations 
is the Fourier transform of the electric fields at aperture, i. e. , z = 0 . 

2.4 Electric Fields in Single Semi- infinite Layer 

When only a single semi- infinite muscle layer is considered, the 
equations in the previous section can be completely utilized by setting 
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z equaL to some convenient value such as one (1) centimeter for 

ijc 5 ' c 

normalization purposes and by Letting = W^, 3 € 2 ‘ 

Another method of obtaining the expressions for a singLe Layer 
medium is to utiLize the equations for region (2) and match boundary 
conditions at z = 0. Both methods are easy to carry out, and since the 
resuLts are reLativeLy simpLe, they wiLL not be incLuded in this discussion. 



2. 5 The Theoreticai Probiem 

An independent choice of any two of the fieid quantities wiLL auto- 

12 

maticaLLy determine aLL other fieid quantities. For purposes of this 
research, the fieLd quantities wiLL be chosen to describe a fieLd configura- 
tion at the aperture which is that corresponding to a TE^ waveguide 

'i 

mode . 

The eLectric fieLd in the x-direction is chosen to be of a uniform 
magnitude and phase along the x-direction and a cosine magnitude varia- 
tion with constant phase in the y-direction. The eLectric fieLd in the 
y-direction is zero everywhere in the aperture. 

This can be mathematicaLLy stated as 



E x (x,y,0) = 0 



-oc < y < 




2. 5. 1 



E (x, y, 0) = Cos — ^ 
x a 



b . . b 

2 ± y - 2 



2 . 5 . 2 



E y ( x > Y> °) 



-oc < x < oo 



0 



2. 5. 3 
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2. 6 Errors Due to Less Than Infinite Integration 

A completely accurate computation of the electric fields at any point 
(x, y, z ) would necessitate integration over the complete infinite u - v 
plane. Such an analysis in this problem is entirely unrealistic and would 
usually make only extremely insignificant contributions outside the normal, 
finite ranges of integration. 

Occasionally, for some integration methods, significant computa- 
tional errors will result when the integration is carried out over finite 
limits. For the problem defined in the previous section, two specific 
instances exist for which finite integration fails to yield answers which 
are devoid of significant computational errors. 

First, when electric fields exist near a conducting wedge (See 
Figure 2. 2), a phenomenon known as edge effect arises. Near a con- 
duction wedge, some field quantities will tend toward infinity as the edge 
is approached even though the energy remains finite. This edge effect 
will produce very high field strengths near the edges of the aperture and 
hot spots in the relative heating pattern. 

1 3 

In the analysis of the edge effect, Guy assumed a single homo- 
geneous semi-infinite lossy medium. This assumption does not account 
for reflections from the fat-muscle interface, and when z , the distance 
from the surface to the fat-muscle interface, is small, this assumption 

introduces another error which must be corrected. 

13 

Guy has analyzed the two problem situations described above and 

derived corrections which are added to E and E when x and z are 

x z 
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in cLose proximity to the aperture edge or when is so small that 
imaging effects must be taken into account. The corrections are given by 




Figure 2 . 2 
The Edge Effect 
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C 
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2 . 6 . 2 



These corrections are to be added to the electric fields at points 
near the edge of the aperture and at ail points in the Region (1) expres- 
sions for E and E when is small enough to cause imaging. Note 

that the integrands in the correction can be subtracted from the integrands 
derived in Section 2. 3 before the integration with respect to u is performed. 

All of the necessary equations for computing relative heating pat- 
terns have now been developed. Attention can now be focused on an 
efficient algorithm to compute the necessary transformations. 



CHAPTER 3 



THE FAST FOURIER TRANSFORM 

3. 1 Background 

Since 1822 when Jean Baptiste Four ier (1768-1830) published his 
Theorie Analytique de la Chaleur , in which he set forth the trigonometric 
series and integrals which are named after him, Fourier analysis has 
played an important role in the analysis of many of the problems of physics 
and engineering. Fourier himself only presented the general boundary 
value problem for the case of heat conduction, but today his methods find 
wide application in the disciplines metnioned above. Besides the heat 
conduction problem where the method is still apropos, many problems of 
today are particularly amendable to the method; to name a few: 

a. Power spectral density 

b. System simulation 

c. Noise analysis 

d. Statistical distribution analysis 

e. Frequency analysis of bandlimited waveforms 

f. Convolution (lagged product) problems 

g. Acoustic wave propogation 

h. Speech transmission 

i. Transport phenomena 

j. Optics 

k. Electro magnetics. 
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3.2 Fourier Analysis 

The theory of Fourier analysis is beyond the scope of this 
14 15 

report, ’ but several excellent texts exist pertaining to this subject. 

The major drawback of Fourier analysis lies in the need to evalu- 
ate integrals which have infinite limits, as is seen in the one - dimensional 
Fourier transform pair utilized in frequency analysis: 



oo 

S(f) I x(t) e " j27Tftdt 
- 00 



3. 2. 1 



x(t) I s(f) ej2Trft df 

-00 



3 . 2.2 



or the equivalent two-dimensional form usually seen in boundary value 
problems.* 

00 

S(u. v) = ~ y ^ d(x, y) e + Vy ^ dxdy 3. 2. 3 

-oo 

oo 

d(x,y) = i J j 1 S(u, v) e j(uX + Vy) dudv 3. 2. 4 

-00 



Equations 3. 2. 3 and 3. 2. 4 are typical of those developed in Chapter 2, 
but for purposes of illustration of the FFT the equations 3.2. 1 and 
3.2.2 will be used. 

3. 3 The Fast Fourier Transform 

Prior to the advent of high speed computers, complex integration 
was handled by such methods as the saddle point or stationary phase 
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techniques which are valid only for expressions of large argument. Thus 
the really important questions of analysis of near zone fields of an antenna 
were seldom investigated since in this area the integrals were unassailable 
by analytic methods, except in a few special cases. The computer changed 
this, but did not significantly advance the state of the art due to the high 
cost of computing time. 



known as the fast Fourier transform. It then became economical to com- 
pute discrete Fourier transforms, and many problems could be economi- 
cally investigated. 



looking at the discrete Fourier transform pair associated with equations 
3. 2. 1 and 3. 2. 2, 



A real advance came when Cooley and Tukey found the algorithm 



Returning now to the Fourier transform itself, and in particular 



N-l 




S(f ) = At x(t ) e 
n k 

k=0 




n = 0, ±1, ±2, ... 



3. 3. 1 




N/2 



k = 0, +1,+ 2 , ... N-l 



3. 3. 2 



n=N/2 



now, making the following substitutions, 



t = kA t 
k 



t = T/N 



f = nAf 



n 




then equations 3. 3. 1 and 3. 3. 2 become, according to Cooley 
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N- 1 



S,„) = £ x(k) <" k »/ N 



k=0 



N-l 



x<k> . £ S,n) J*' ( " k) / N 



n=0 



3. 3. 3 



n = 0, + 1 , +2, . . . N-l 
k = 0, + 1,42, ... N-l 



3. 3. 4 



3. 4 Factoring the Transform 

For purposes of illustration, when equation 3. 3. 3 is written out 

i 

in full for N = 4 and with W = e 
define 



l -j2irf/N . , 18 

r - ° , Brigham and Morrow 



[S(n)j = 



[ S(n)J = 



r ° 


0 


w 


W 1 


, 0 


1 


w 


W 1 


0 


2 


w 


W 1 


0 


3 


w 


W 1 


n)] as: 


S(0) 




S(2) 




S(1) 




S(3_)_ 





0 0 
r W 

2 3 

r W 

4 6 

r W 

r 6 „ r 9 



x 0 (°) 

x o (1) 

x o (2) 

x 0 (3) 



0 0 0 0 

w w w w 

0 2 4 6 

W W W W 

0 12 3 

W W W W 

0 3 6 9 

W W W W 



S(0) 

S(l) 

S(2) 

S(3) 



3. 4. 1 



x 0 (°) 

x o (1) 

x o (2) 

_ x o (3 2j 



= [w] [x 0 (n)] 



3. 4. 2 



r u 

The [ W ] matrix in 3. 4. 2 can now be factored, and since W =1 
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[ S(n)] = 



1 W 0 0 



1 W 0 0 



0 0 1 W 



0 0 1 W' 



1 1 W° 0 



0 10 w 



low 0 



0 10 w 



Alternately, this can be written: 



1 W 0 0 



1 W 0 0 



[ S(n) ] = [x (n)] = 

2 



0 0 1 w 



0 0 1 W' 



where: 



[^(n)] = 



10W 0 



0 10 W 



1 0 W 0 



0 10 W. 



x 0 (°) 

x o (1) 

x o (2) 

l x o (3 2i 



[ X x ( n ) ] 



t x 0 ( n )] 



3. 4. 3 



3. 4. 4 



3- 4. 5 



This factorization looks somewhat mystical, so to put minds at 
ease and insure that the factorization is valid, and to insure that equation 



3. 4. 4 agrees with 3. 3. 3, from equation 3. 4. 4, 

Xj(o) + x^nw 0 " 1 



Xl (o) + x^nw 

X x (2) + Xl (3)W ] 
Xj(2) + Xl (3)W : 



/ 



[ S(n) ] = 



3. 4. 6 



and from equation 3. 4. 5, 
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[x^n)] 



x Q (0) + x Q (2)W 
x 0 (l) + x Q (3)W ( 
x Q (0) + x Q (2)w‘ 

,x 0 (l) + x Q (3)W‘ 



3. 4. 7 



then, substituting 3.4.7 into 3.4.6 yields 



[ S( n) ] = 



X 1 

O 

O 


+ 


0 

x o (2)W 


+ 


< x o (1) 


+ 


0 0“| 

x Q (3)W )W 


x 0 (°) 


+ 


x 0 (2)W° 


+ 


(x 0 (D 


+ 


x 0 (3)W°)W 2 


x 0 (°) 


+ 


x Q (2)W 2 


+ 


<*o< l > 


+ 


x q (3)W 2 )W 1 


- x o< 0) 


+ 


x Q (2)W 2 


+ 


< x 0* l) 


+ 


x q (3)W 2 )W 3 _ 



3. 4. 8 



n n±4 n n±4k 

If we note that W = W or, more generally speaking W = W , 



where k is some integer, then 3. 4. 8 upon rearrangement becomes 



once again equation 3. 4. 2 written out in full, instead of matrix form: 



[ S(n) ] = 



x Q (0) + x 0 (l)W° + x 0 (2)W° + x 0 (3)W° 

x 0 (0) + x 0 (l)W 2 + x Q (2)W 4 + x Q (3)W 6 

x Q (0) + x Q (l)W 1 + x Q (2)W 2 + x Q (3)W 3 

x Q (0) + x Q (l)W 3 + x Q (2)W 6 + x Q (3)W^ 



3. 4. 9 



One discrepancy has arisen in that we have derived [ S(n) ] instead 
of S(n). This can very easily be corrected by an algorithm which looks 
at n as a binary address of a particular element of the matrix S. To 
find the correct address is simply a matter of bit reversing the binary 



number n. 
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Then, 



[ S(n)] = 



S(00) 

S(01) 

S(10) 

S(ll) 



and bit reversing the address yields 



[ S(n)] 



S(00) 

S(10) 

S(01) 

S(ll) 



3. 4. 10 



3. 4. 11 



While it is true that the above development applies only to the 
case N = 4 and not to the general case, when the discrete Fourier 
transform, equation 3. 3. 3, is examined for a general case, as has been 

19 

shown by Gentleman and Sande, then the factorization can be performed 
for any N that is not a prime number. This factorization is a recursive 
formula which defines larger Fourier transforms in terms of smaller 
transforms. A multi-dimensional discrete Fourier transform, such as 
this research investigates, can be efficiently computed by performing 
the factorization in each dimension separately. 



3. 5 The Algorithm 

The obvious question to ask now, is how to obtain the factored 
matrices indicated in equation 3. 4. 3. This is most easily seen in a 
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tree graph (See Figures 3. 1 and 3. 2) where a vertical column of nodes 

will correspond to the x^(k) arrays (where for our example & = 0, 1,2). 

For cases where N = 2^(y an integer), we will have y + 1 vertical arrays 

2 

and N horizontal arrays. For this example we will use y = 2, N ^ 2 = 4. 

In the tree graph, a solid line will indicate that the node in the previous 
vertical array is to be multiplied by the W factor raised to the power 
indicated in the node of interest, then added to the value of the node from 
which the dashed line originates. The following algorithm yields the neces- 
sary node integers. 

1. Express the argument k as a binary number. 

th 

2. For the S. array, shift the number y - i places to the 
right and reverse the bits. 

To find out where the solid and dashed lines originate: 

1. Let the argument k be expressed as a binary number such 



that k - k _ , k _, . . . k_ . Thus if k = 6, y = 3 then 
y- 1 y-2 0 

k 2 = 1. k ! = 1. k Q = 0 or k 2 k i k () = H0) 2 = 6)^. 

th th 

2. Now in the i array, node k has a solid line from k 

node in the S. - 1 array except that k ^ is a one (1). 

3. A dashed line originates from address k except that 

k . is a zero (0). 

Y--* 

To see how this works, apply the rules to two cases N 4, N = 8 
and arrive at the tree graphs shown in Figures 3. 1 and 3. 2 respectively. 
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K L= I L= 2 

n 




.►Complex Addition 
.►Complex Multiplication 



N=4 y-2 



FIGURE 3. 1 



Tree Graph, N = 4 
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N = 8 S 2 3 



/=3 



K 

ii 



L s 



L= 2 



L= 3 




.Complex Addition 
^.Complex Multiplication 



FIGURE 3. 2 



Tree Graph, N = 8 
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Now verify the tree graph by looking at S(l) = x (1) in Figure 

c* 



3. 1. 



x 2 (l) = Xl (2) + x^JW 1 
x^ (2) = x 0 (0) + x Q (2)W 2 
x x (3) = x 0 (l) + x Q (3)W 2 



3. 5. 1 



Then, 



S(1 ) = x 2 (l) = x Q (0) + x Q (2)W 2 + (x 0 (l) + x () (3)W 2 )W 1 
S( 1 ) = x Q (0) + x Q (l)W 1 + x Q (2)W 2 + x q (3)W 3 , 



3. 5. 2 



which is in complete agreement with prior results. 

3. 6 Advantage of the FFT 

The conventional method of performing Fourier analysis on a set 
of data points consists of summing a finite number of terms in a Fourier 
series to obtain each transformed data point. Each term added into the 
series requires a complex multiplication to obtain the term and a com- 
plex addition to include the term in the series. Thus, if N data points 
are to be transformed utilizing N terms of the series and one complex 

multiplication and one complex addition constitute one operation, then 

2 

N operations must be performed to accomplish the transformation. 

4 

Likewise, a two-dimensional array of data points will require N 
operations. 



By comparison, the FFT requires N operations for each set 
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Nlog^N operations for a one-dimensional transform* For the 

2 

two-dimensional case 4N log^N operations are required. 

With the above results, the time savings factor (TSF) for a 
two-dimensional array of 64 x 64 points is 



TSF = 



4N log N 41og N 

Cd Cd x 



N' 



N 



170 



3. 6. 1 



3. 7 Problems in FFT Analysis 

Aliasing, leakage and the "picket fence" effect are the most com- 

20 

mon errors encountered in utilization of the fast Fourier transform. 

The effects of these errors can usually be completely eradicated, or at 
the very least be reduced to acceptable levels, by insuring that the data 
is sampled at least twice during each cycle and by proper selection of 
data window. 



3. 8 Space and Frequency Relationships 

A problem encountered in performing either the Fourier trans- 
form or the inverse Fourier transform is the correlation of a sampling 
rate in the frequency domain to data point incrementation in the time or 
special domain. In conventional Fourier analysis, the problem is 
resolved by the very nature of the series, in that the argument of the 
functions in each term of the Fourier series contains time and frequency 
explicitly as seen in equation 3. 8. 1. 



33 



oc 



“> * Z c n 





3. 8. 1 



The Fourier series does not demand the computation of x at any particu- 
lar t. The onLy requirement is that t lie within the data window 
-T/2 < t < T/2 in the asymmetric case and within 0 < t < T in the 
symmetric case when only N terms of the series are utilized. 

By contrast, time and frequency are nowhere to be found in the 
formulation of the fast Fourier transform. This can lead to confusion 
and uncertainty regarding the exact coordinates of a point in either domain 
unless careful attention is paid to equation 3. 8. 2. 



These relations are simple, but can mislead the user of the FFT. 

In order to clear up the question surrounding this problem, con- 
sider a one -dimensional problem in which an array of data points are to 
be (1) transformed into the frequency domain; (2) multiplied by some 
other function; (3) then the product transformed back into the time domain. 
The problem may be stated as 



1 



2tt 



3. 8. 2 



1 NAf 



NAu 



N-l 




S(k) = 



3. 8. 3 
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N-l 



x i< n) ' TJ 



S(k)G(k) e' 



j 2tt nk/N 



k=0 



3. 8. 4 



with the following definitions (See Figure 3. 3). 

T = width of data window 
N = number od data points 

These independent choices dictate the following relations. 



Ax 



0 



T 

N-l 



Aou 



2 it (N-l) 
NT 



F 



= Spectral window 



N-l 

NT 



Ax 



l 



T 

N-l 



3. 8. 5 



3. 8. 6 



Then, assuming G is a continuous function, it must be evaluated 
in the interval -F/2=w = F/2 at N points spaced equally in increments 
of A o>. This also assumes that G is a function of o> rather than f, 

and this is quite often the case, e. g. , the relative heating equations in the 
previous chapter which have been programmed utilizing the FFT and 
the results compared with a Gaussian quadrature transformation in the 



next chapter. 



CHAPTER 4 



UTILITZATION OF THE FFT ON RELATIVE HEATING PROBLEMS 



4. 1 Introduction 



Diathermy, the therapeutic heating of deep body tissue, is an 
important part of physical medicine and rehabilitation. Research at the 
University of Washington in this field is concerned with the determina- 
tion of an optimum frequency and aperture configuration which will 
yield good heating in musculature without excess heating in the sub- 
cutaneous fat. This can be recognized in a relative heating pattern 

\ 

which shows low heating in the fat relative to the amount of heating 
observed in the muscle. 

The research basically has two phases. The first stage is a 
theoretical analysis of the problem utilizing plastic models of tissue 
and computer models of the electromagnetics. The second phase is the 
building and testing of diathermy applicators, once again using the 
already proven plastic and computer models. 

Chapters 2 and 3 have laid the groundwork for a computer model 
of the diathermy problem. The succeeding sections of this chapter 
will finish the construction of the computer model. Each phase of the 
research has a similar, but different computer model and each will be 



discussed separately. 
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4.2 The Theoretical Phase 

The first phase computer model for diathermy research consists 
of an executable computer algorithm which can evaluate the relative 
heating expressions of Chapter 2 for a theoretical aperture distribution. 
A suitable algorithm which can accomplish the transformations is the 
fast Fourier transform. 

Several versions of the fast Fourier transform can be found in 

21 , 22 

the literature, some of which are written in Fortran and others 

21 

are written in Algol. The algorithm of Eguchi has been used in this 
research for two reasons; 

a. ) its language is Fortran, 

b. ) it can be implemented on the University of Washington 

computers without extensive revision. 

Some revision of Eguchi’ s program has been made to cut out the 
portions of his program which are not necessary in the diathermy prob- 
lem. Since extensive calculations must be done to build the arrays to 
be transformed and the amount of memory used is a determining factor 
in computer costs, features in Eguchi' s program such as magnitude 
distributions, convolutions, gradients and the various output options 
are unusable and have been discarded since they obligate memory that 
can be more efficiently used. 

Figure 4, 1 gives a flow chart for the use of the various sub- 
routines given in the Appendices and outlines the method of computer 
solutions necessary in the first phase of diathermy research. 
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FIGURE 4. 1 



Flow Chart of Theoretical Problem 
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h 


1 













00 81 
INDEX =1,3 



V= VINIT 
DO 45 
M = 1,64 



U= UINIT 
DO 44 





FIGURE 4, 1 (continued) 
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FIGURE 4. 1 (continued) 
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FIGURE 4. 1 (continued) 
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4. 3 The Empirical Phase 

Upon successful determination of an optimum frequency and 
aperture configuration, extensive tests must be made of any prototype 
diathermy applicators. 

Utilizing the methods of the first phase computer solution and 
having available a sufficient amount of empirical data from the field 
configuation at the surface of the plastic models, the fast Fourier 
transform can once again be used to determine the relative heating 
patterns throughout the volume of the planar model. 

A suitable method of solution for problems of an empirical 
nature is outlined in Figure 4. 2. 

4.4 Objective Restatement 

This research is primarily concerned with evaluating the relative 
merits of the fast Fourier transform to the solution of the theoretical 
relative heating problem versus the Gaussian quadrature method 
currently in use. Solutions for empirical problems will follow as a 
direct by-product of this research. The results of computations by 
both methods are compared in the next chapter. 
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FIGURE 4. 2 



Flow Chart of Empirical Problem 
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FIGURE 4. 2 (continued) 
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FIGURE 4. 2 (continued) 
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FIGURE 4. 2 (continued) 



CHA PTER 5 



RESULTS 

5. 1 Introduction 

Relative heating patterns (RHPS) as computed by the fast Fourier 
transform (FFT) and Gaussian quadrature (GQ) methods will be presented 
in the succeeding sections of this chapter. An error criterion will be 
defined and the computed relative heating values will be compared on a 
point by point basis whenever possible. 

5. 2 The Error Criterion 

Any digital integration method has an associated, inherent error. 

Quite often, the exact error cannot be ascertained; i. e. , the bounds of the 

error are all that can be established. At other times, even the error 

bounds are difficult to establish. In the case of the FFT, Gentleman and 
13 

Sande have proven the bounds of error due to round off, but nowhere is 
there a complete proof of error bounds for a fast Fourier transformation. 
No attempt will be made here to set error bounds. Instead, the FFT will 
be compared to the best available values from the Gaussian quadrature 
method. 

Since both methods have inherent inaccuracies due to less than 
infinite integration, any error encountered in comparing results will also 
contain inaccuracies which are impossible to establish. For purposes of 
this research, error will be defined as the absolute difference in relative 
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heating values at the point (x,y, z) as computed by the Gaussian quadrature 
and fast Fourier transform methods. This error will be put on a percentage 
basis by dividing the magnitude of the error by the magnitude of the nor- 
malization point always has a value of unity in both methods, the point 
(0, 0, z^) will always have zero error when comparing the two transforma- 
tion methods. This definition is then 

% ERROR = 100 (PGAUSS - PFFT) 5.2. 1 

where PGAUSS and PFFT are the relative heating values computed at 
(x,y,z) by the Gaussian quadrature and fast Fourier transformations 
r espectively. 

To be an acceptable tool for computing relative heating patterns, 
the fast Fourier transform must compute values to within 5% ERROR, 
except at points which lie along lines with very steep slopes. 

5. 3 Data Presentation 

Relative heating values and the associated error are tabulated in 
Tables 5.2 - 5.6 for several different frequencies. Each table is fol- 
lowed by CALCOMP plots which give a visual representation of the rela- 
tive heating patterns as computed by the FFT and Gaussian quadrature 
methods, in that order. Table 5. 1 list the frequencies, associated 
dielectric constants, and loss tangents used in each of the five cases 
considered. East case considers a rectangular aperture (12 x 16cm. ) 
in an infinite ground plane (See Figure 2. 1). The relative heating 
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patterns were computed at Z = 0. 1, 1. 0, 2. 0, 3. 0, 4. 0cm. with the 

fat muscle interface occurring in each case at = 2. 0cm. 

A brief explanation is in order concerning the CALCOMP plots 
which follow each tabulated comparison of the two methods. The first 
plot is a relative heating pattern as computed by the FFT. This is 
followed immediately by a transparency taken from the FFT results. 

This transparency overlays a relative heating pattern for the same case 
as computed by the GQ method. 

The transparencies unavoidably shrunk somewhat in processing, 
nevertheless they provide an excellent visual comparison of the two com- 
putational methods. For example, examine Figures 5.3, 5.3a, 5.4 
which are the relative heating patterns for the 7 50 mHz comparison in 
Table 5. 3. The following agreements and discrepancies are readily 
identifiable in the overlay: 

a. There is very good agreement in the two plots for Z = 1 
and 4 cm. 

b. The Z = 3 cm. curves differ by about 1% in the range 
3 cm. < x < 8 cm. 

c. About 2% difference can be seen in the Z = 2 cm. curves 
in the range 1 cm. < x < 8 cm. 

d. The FFT curve for Z = 0. 1 cm. has some aliasing due to the 
truncated series, but shows good agreement with the GQ curve 
except at the edge of the aperture. 

Similar comparisons can be made for each frequency, but will not be made 



here. 
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Table 5. 1 
Input Data 



F 




6 

m 


T 

m 


e f 


T f 


400 


MHZ 


52. 79 


0. 9842 


5. 55 


0. 3745 


750 


MHZ 


51. 54 


. 5944 


5. 59 


. 2659 


918. 8 


MHZ 


51. 4 


. 488 


5.75 


. 25 


1500 


MHZ 


49- 52 


. 3858 


5. 58 


. 1943 


2450 


MHZ 


47. 42 


0. 3059 


5. 56 


0. 1541 




F = 


frequency 










e = 

m 


dielectric constant 


in muscle 








T = 
m 


loss tangent in muscle 








e r = 


dielectric constant 


in fat 







= loss tangent in fat 
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Table 5.2 
400 MHZ 



Comparison of RHPS as Computed by FFT and GQ 



Z = 0. 1 cm Y = 0. 0 cm Z ^ =2.0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


0. 0 


0. 169388 


0. 169952 


0. 056 


0. 5 


. 175147 


. 170537 


. 461 


1. 0 


. 171048 


. 172285 


. 124 


1. 5 


. 178611 


. 175194 


. 342 


2. 0 


. 175963 


. 179276 


. 331 


2. 5 


. 185448 


. 184532 


. 092 


3. 0 


. 183975 


. 190959 


. 698 


3. 5 


. 195594 


. 198569 


0. 297 


4. 0 


. 194992 


. 207351 


1. 236 


4. 5 


. 209297 


.217296 


0. 800 


5. 0 


. 209303 


. 228461 


1. 916 


5. 5 


. 227837 


. 241038 


1. 320 


6. 0 


. 228445 


.255586 


2.714 


6. 5 


.256536 


. 273632 


1.710 


7. 0 


. 259753 


. 299380 


3. 963 


7. 5 


. 332851 


. 347174 


1. 432 


8. 0 


. 326463 


. 437956 


11. 149 


8. 5 


. 151589 


. 190027 


3. 844 


9- 0 


. 090918 


. 130376 


3. 946 


9. 5 


. 066715 


. 100687 


3. 397 


10. 0 


0. 051848 


0. 081288 


2. 944 




Z = 1. 0 cm Y = 0. 


0 cm Z^ =2.0 cm 




0. 0 


0. 091048 


0. 090634 


0. 041 


0. 5 


. 091409 


. 091156 


. 025 


1. 0 


. 092435 


. 092712 


. 028 


1. 5 


. 094167 


. 095286 


. 112 


2. 0 


. 096532 


. 098847 


. 231 


2. 5 


. 099564 


. 103357 


. 379 


3. 0 


. 103185 


. 108777 


. 559 


3. 5 


. 107442 


. 115069 


. 763 


4. 0 


. 112266 


. 122200 


0. 993 


4. 5 


. 117727 


. 130148 


1. 242 


5. 0 


. 123753 


. 138889 


1. 514 


5. 5 


. 130372 


. 148365 


1.799 


6. 0 


0. 137287 


0. 158386 


2. 110 
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Table 5. 2 (continued) 





Z = 1. 0 cm Y = 


0. 0 


cm Z^ =2.0 cm 




X 


P FFT 




P GAUSS 


% ERROR 


6. 5 


0. 143945 




0. 168356 


2. 441 


7. 0 


. 148450 




. 176575 


2. 812 


7. 5 


. 146659 




. 178684 


3. 203 


8. 0 


. 132127 




. 167522 


3. 539 


8. 5 


. 107165 




. 143679 


3. 651 


9. 0 


. 083202 




. 118478 


3. 528 


9. 5 


. 065198 




. 097531 


3. 233 


10. 0 


0. 052651 




0. 080851 


2. 820 




Z = 2. 0 cm Y = 


0. 0 


cm Zj =2.0 cm 




0. 0 


1. 000000 




1. 000000 


0. 000 


0. 5 


0. 997064 




0. 996597 


. 047 


1. 0 


. 988276 




. 986425 


. 185 


1. 5 


. 973692 




. 969587 


. 41 1 


2. 0 


. 953406 




. 946253 


0. 715 


2. 5 


. 927548 




. 916656 


1. 089 


3. 0 


. 896284 




. 881081 


1. 520 


3. 5 


. 859815 




. 839869 


1. 995 


4. 0 


.818375 




.793416 


2. 496 


4. 5 


. 772241 




. 742172 


3. 007 


5. 0 


. 721731 




. 686660 


3. 507 


5. 5 


. 667224 




. 627489 


3. 974 


6. 0 


. 609190 




. 565400 


4. 379 


6. 5 


. 548269 




. 501325 


4. 694 


7. 0 


. 485449 




. 436549 


4. 890 


7. 5 


. 422356 




. 372947 


4. 941 


8. 0 


. 361438 




. 313075 


4. 836 


8. 5 


. 305438 




.259591 


4. 585 


9. 0 


. 256178 




.214073 


4. 210 


9- 5 


. 213941 




. 176540 


3. 740 


10. 0 


0. 177985 




0. 146014 


3. 197 




Z = 3. 0 cm Y = 


0. 0 


cm Z i =2.0 cm 




0. 0 


0. 550351 




0. 551602 


0. 125 


0. 5 


. 548625 




. 549611 


. 099 


1. 0 


. 543457 




. 543656 


. 020 


1. 5 


. 534874 




. 533789 


. 109 


2. 0 


. 522926 




. 520097 


. 283 


2. 5 


. 507689 




. 502710 


. 498 


3. 0 


0. 489268 




0. 481801 


0. 747 
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Table 5. 2 (continued) 



Z = 3. 0 cm Y - 0. 0 cm Z ^ =2.0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


3. 5 


0. 467811 


0. 457597 


1. 022 


4. 0 


. 443515 


. 430385 


1. 313 


4. 5 


. 416643 


. 400537 


1.611 


5. 0 


. 387545 


. 368519 


1. 903 


5. 5 


. 356670 


. 334904 


2. 177 


6. 0 


. 324587 


. 300397 


2. 419 


6. 5 


. 291979 


. 265826 


2. 615 


7. 0 


. 259628 


. 232110 


2. 752 


7. 5 


. 228338 


. 200174 


2. 816 


8. 0 


. 198829 


. 170832 


2. 800 


8. 5 


. 171621 


. 144637 


2. 698 


9- 0 


. 146975 


. 121814 


2. 516 


9- 5 


. 124916 


. 102294 


2. 262 


10. 0 


0. 105311 


0. 0858 08 


1. 950 




Z = 4. 0 cm Y = 0. 


0 cm Z ^ =2.0 cm 




0. 0 


0. 300086 


0. 301578 


0. 149 


0. 5 


. 299063 


. 300403 


. 134 


1. 0 


. 296001 


. 296889 


. 089 


1. 5 


. 290926 


. 291075 


. 015 


2. 0 


. 283882 


. 283025 


. 086 


2. 5 


. 274937 


. 272840 


. 210 


3. 0 


. 264192 


. 260654 


. 3 54 


3. 5 


. 251782 


. 246649 


. 513 


4. 0 


. 237889 


. 231057 


. 683 


4. 5 


. 222742 


.214170 


0. 857 


5. 0 


. 206626 


. 196335 


1. 029 


5. 5 


. 189871 


. 177952 


1. 192 


6. 0 


. 172844 


. 159463 


1. 338 


6. 5 


. 155923 


. 141324 


1. 460 


7. 0 


. 139465 


. 123966 


1. 550 


7. 5 


. 123766 


. 107756 


1. 601 


8. 0 


. 109037 


. 092960 


1. 608 


8. 5 


. 095389 


. 079721 


1. 567 


9. 0 


. 082849 


. 068066 


1. 478 


9- 5 


. 071386 


. 057938 


1. 345 


10. 0 


0. 060946 


0. 049218 


1. 173 



RELRTIVE HEATING IN X-Z PLANE, 




400 mHz FFT RHP 



RELATIVE HEATING IN X-Z PLANE, 




400 mHz GQ RHP 
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Table 5. 3 
750 MHZ 



Comparison of RHPS as Computed by FFT and GQ 



Z = 0. 1 cm Y = 0. 0 cm Z^ =2.0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


0. 0 


0. 219999 


0. 230121 


1. 012 


0. 5 


. 227923 


. 230939 


0. 302 


1. o 


. 223827 


. 233344 


. 952 


1. 5 


. 235375 


. 237112 


. 174 


2. 0 


. 233776 


. 241917 


.814 


2. 5 


. 247462 


. 247393 


. 007 


3. 0 


. 246062 


.253100 


. 704 


3. 5 


. 260112 


.258668 


. 144 


4. 0 


. 256719 


. 263765 


. 705 


4. 5 


. 270403 


. 268169 


. 223 


5. 0 


. 263880 


. 271879 


. 800 


5. 5 


. 278808 


. 275113 


. 370 


6. 0 


. 269301 


. 278575 


. 927 


6. 5 


. 291805 


. 283724 


0. 808 


7. 0 


. 280616 


. 293975 


1. 336 


7. 5 


. 345529 


. 320156 


2. 537 


8. 0 


. 298479 


. 353216 


5. 474 


8. 5 


.214655 


. 144184 


1. 047 


9- 0 


. 073904 


. 071551 


0. 235 


9- 5 


. 052573 


. 052983 


0. 041 


10. 0 


0. 037884 


0. 041624 


0. 374 




Z = 1 . 0 cm Y = 0. 


0 cm Zj =2.0 cm 




0. 0 


0. 111816 


0. 113543 


0. 173 


0. 5 


. 112691 


. 114293 


. 160 


1. 0 


.115150 


. 116457 


. 131 


1. 5 


. 118966 


. 119783 


. 082 


2. 0 


.213633 


. 123896 


. 026 


2. 5 


. 128689 


. 128342 


. 035 


3. 0 


. 133528 


. 132638 


. 088 


3. 5 


. 137741 


. 136397 


. 134 


4. 0 


. 140919 


. 139277 


. 164 


4. 5 


. 142974 


. 141132 


. 184 


5. 0 


. 143856 


. 141959 


. 190 


5. 5 


. 143801 


. 141860 


. 194 


6. 0 


. 142846 


. 140888 


. 196 


6. 5 


0. 140843 


0. 138709 


. 213 
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TabLe 5. 3 (continued) 



Z = 1. 0 cm Y = 0. 0 cm Z ^ = 2. 0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


7. 0 


0. 136485 


0. 133994 


0. 249 


7. 5 


. 216978 


. 123763 


. 321 


8. 0 


. 1091 12 


. 105312 


. 380 


8. 5 


. 086657 


. 083343 


. 331 


9. 0 


. 066680 


. 065268 


. 141 


9. 5 


. 050719 


. 051936 


. 122 


10. 0 


0. 037949 


0. 041896 


0. 395 




Z = 2. 0 cm Y = 0. 0 


cm Z^ + 2. 0 cm 




0. 0 


1. 000000 


1. 000000 


0. 000 


0. 5 


0. 993348 


0. 993783 


. 043 


1. 0 


. 973627 


. 975335 


. 171 


1. 5 


. 941525 


. 945266 


. 374 


2. 0 


. 898147 


. 904547 


. 640 


2. 5 


. 844959 


. 854473 


0. 951 


3. 0 


. 783711 


. 796589 


1. 288 


3. 5 


. 0716354 


. 732627 


1. 627 


4. 0 


. 644953 


. 664420 


1. 947 


4. 5 


. 571594 


. 593836 


2. 224 


5. 0 


. 498302 


. 522695 


2. 439 


5. 5 


. 426970 


. 452710 


2. 574 


6. 0 


. 359307 


. 385434 


2. 613 


6. 5 


.296831 


. 322244 


2. 541 


7. 0 


. 240923 


. 264404 


2. 348 


7. 5 


. 192899 


. 213171 


2. 027 


8. 0 


. 153845 


. 169734 


1. 589 


8. 5 


. 124032 


. 134680 


1. 065 


9. 0 


. 102427 


. 107445 


0. 502 


9. 5 


. 087076 


. 086560 


0. 052 


10. 0 


0. 075927 


0. 07 0414 


0. 551 




Z = 3. 0 cm Y = 0. 0 


cm Z^ =2.0 cm 




0. 0 


0. 520233 


0. 519987 


0. 025 


0. 5 


. 516642 


. 516633 


. 001 


1. 0 


. 506000 


. 506682 


. 068 


1. 5 


. 488694 


. 490472 


. 178 


2. 0 


. 465339 


. 468543 


. 320 


2. 5 


. 436742 


. 441603 


. 486 


3. 0 


. 403861 


. 410496 


. 663 


3. 5 


. 367755 


. 376162 


0. 841 


4. 0 


. 329544 


. 339601 


1. 006 


4. 5 


0. 290378 


0. 301852 


1. 147 
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Table 5. 3 (continued) 



Z = 3. 0 cm Y = 0. 0 cm - 2. 0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


5. 0 


0. 251410 


0. 263964 


1. 255 


5. 5 


. 213774 


. 226984 


1. 321 


6. 0 


. 178563 


. 191927 


1. 336 


6. 5 


. 146774 


. 159773 


1. 296 


7. 0 


.119219 


. 131172 


1. 195 


7. 5 


. 096394 


. 106731 


1. 034 


8. 0 


. 078352 


. 086502 


0. 815 


8. 5 


. 064673 


. 070173 


. 550 


9- 0 


. 054853 


. 057157 


.2 57 


9. 5 


. 047166 


. 046785 


. 038 


10. 0 


0. 041558 


0. 038461 


0. 310 




Z = 4. 0 cm Y = 0. 


0 cm =2.0 cm 




0. 0 


0. 270218 


0. 270042 


0. 018 


0. 5 


. 268288 


. 268239 


. 005 


1. 0 


. 262569 


. 262887 


. 032 


1. 5 


. 253261 


. 254162 


. 090 


2. 0 


. 240690 


. 242347 


. 166 


2. 5 


. 225288 


. 227822 


. 253 


3. 0 


.207581 


.211 048 


. 347 


3. 5 


. 188166 


. 192560 


. 439 


4. 0 


. 167699 


. 172950 


. 525 


4. 5 


. 146866 


. 152844 


. 598 


5. 0 


. 126363 


. 132888 


. 652 


5. 5 


. 106861 


.113707 


. 685 


6. 0 


. 088956 


. 095864 


. 691 


6. 5 


. 073124 


. 079805 


. 668 


7. 0 


. 059664 


. 065816 


. 615 


7. 5 


. 048675 


. 053989 


. 531 


8. 0 


. 040046 


. 044232 


. 419 


8. 5 


. 033496 


. 036313 


. 282 


9. 0 


. 028636 


. 029927 


. 129 


9. 5 


. 025039 


. 024768 


. 027 


10. 0 


0. 022295 


0. 020567 


0. 173 



RELATIVE HEATING IN X-Z PLANE , 




o 



o o o 



750 mHz FFT RHP 



RELATIVE HEATING IN X-Z PLANE - 



o 

II 

>- 
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o 



o 



o 



o 



750 mHz GQ RHP 



60 



Table 5.4 
918. 8 MHZ 

Comparison of RHPS as Computed by FFT and GQ 



Z = 0. 1 cm Y = 0. 0 cm Z ^ = 2. 0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


0. 0 


0. 307665 


0. 319533 


1. 187 


0. 5 


. 318697 


. 320667 


0. 197 


1. o 


. 312688 


. 323887 


1. 120 


1. 5 


. 328181 


. 328736 


0. 055 


2. 0 


. 324581 


. 334546 


. 997 


2. 5 


. 341489 


. 340525 


. 096 


3. 0 


. 336246 


. 346964 


. 972 


3. 5 


. 351737 


. 350319 


0. 142 


4. 0 


. 342179 


. 353316 


1.114 


4. 5 


. 356449 


. 355085 


0. 136 


5. 0 


. 342905 


. 356094 


1. 319 


5. 5 


. 360407 


. 357217 


0. 319 


6. 0 


. 345085 


. 359791 


1. 471 


6. 5 


. 375579 


. 365824 


0. 976 


7. 0 


. 359559 


. 379465 


1. 991 


7. 5 


. 445967 


. 413545 


3. 242 


8. 0 


. 369174 


. 447123 


7. 795 


8. 5 


. 147920 


. 132507 


1. 541 


9. 0 


. 092290 


. 080835 


1. 146 


9. 5 


. 070209 


. 059541 


1. 067 


10. 0 


0. 054302 


0. 046853 


0. 745 




Z = 1. 0 cm Y = 


0. 0 cm Z i =2.0 cm 




0. 0 


0. 144381 


0. 146859 


0. 248 


0. 5 


. 145761 


. 148084 


. 232 


1. 0 


. 149585 


. 151553 


. 197 


1. 5 


. 155298 


. 156693 


. 139 


2. 0 


. 161878 


. 162667 


. 079 


2. 5 


. 168374 


. 168539 


. 017 


3. 0 


. 173745 


. 173442 


. 030 


3. 5 


. 177413 


. 176752 


. 066 


4, 0 


. 179011 


. 178196 


. 082 


4. 5 


. 178812 


. 177896 


. 092 


5. 0 


. 177246 


. 176304 


. 094 


5. 5 


. 175111 


. 174003 


.111 


6. 0 


. 172744 


. 171344 


. 140 


6. 5 


. 169904 


. 167855 


. 205 


7. 0 


. 164432 


. 161370 


. 306 


7. 5 


. 152263 


. 147447 


. 482 


8. 0 


0. 130205 


0. 122927 


0. 728 
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Table 5.4 (continued) 



Z = 1. 0 cm Y = 0. 0 cm Z ^ =2.0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


8. 5 


0. 105004 


0. 095295 


0. 971 


9- 0 


. 084526 


. 073889 


1. 064 


9. 5 


. 068273 


. 058699 


0. 957 


10. 0 


0. 054134 


0. 047413 


0. 672 




Z -2.0 cm Y = 


0. 0 cm Z^ =2.0 cm 




0. 0 


1. 000000 


1. 000000 


0. 000 


0. 5 


0. 994130 


0. 994035 


. 010 


1. 0 


. 976754 


. 976381 


. 037 


1. 5 


. 948539 


. 947722 


. 082 


2. 0 


. 910485 


. 909097 


. 139 


2. 5 


. 863800 


. 861776 


. 202 


3. 0 


. 809766 


. 807131 


. 264 


3. 5 


. 749643 


. 746547 


. 310 


4. 0 


. 684639 


. 681373 


. 327 


4. 5 


.615931 


. 612953 


. 298 


5. 0 


. 544743 


. 542659 


. 208 


5. 5 


. 472431 


. 471968 


. 046 


6. 0 


. 400548 


. 402477 


. 193 


6. 5 


. 330890 


. 335891 


. 500 


7. 0 


. 265553 


. 274031 


0. 848 


7. 5 


. 206994 


.218845 


1. 185 


8. 0 


. 157785 


. 172206 


1. 442 


8. 5 


. 119638 


. 135132 


1. 549 


9. 0 


. 092391 


. 107008 


1. 462 


9- 5 


. 074256 


. 085977 


1. 172 


10. 0 


0. 062899 


0. 070018 


0. 712 




Z = 3. 0 cm Y = 


0. 0 cm Z^ = 2. 0 cm 




0. 0 


0. 521890 


0. 520860 


0. 103 


0. 5 


. 518547 


. 517491 


. 106 


1. o 


. 508664 


. 507526 


. 114 


1. 5 


. 492652 


. 491379 


. 127 


2. 0 


.471135 


. 469684 


. 145 


2. 5 


. 444872 


. 443222 


. 165 


3. 0 


. 414672 


. 412841 


. 183 


3. 5 


. 381339 


. 379406 


. 193 


4. 0 


. 345638 


. 343759 


. 188 


4. 5 


0. 308315 


0. 306731 


0. 158 
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Table 5. 4 (continued) 

\ 



Z = 3. 0 cm Y - 0. 0 cm - 2» 0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


5. 0 


0. 270145 


0. 269174 


0. 097 


5. 5 


. 231993 


. 232008 


. 001 


6. 0 


. 194862 


. 196241 


. 138 


6. 5 


. 159885 


. 162945 


. 306 


7. 0 


. 128228 


. 133128 


. 490 


7. 5 


. 100900 


. 109544 


. 664 


8. 0 


. 078529 


. 086485 


.796 


8. 5 


. 061217 


. 0697 07 


. 849 


9- 0 


. 048567 


. 056560 


.799 


9- 5 


. 039859 


. 046257 


. 640 


10. 0 


0. 034219 


0. 038083 


0. 386 




Z = 4. 0 cm Y = 0. 


0 cm Z i =2.0 cm 




0. 0 


0. 272516 


0. 271596 


0. 092 


0. 5 


. 270639 


. 269714 


. 093 


1. 0 


. 265092 


. 264150 


. 094 


1. 5 


. 256111 


.255142 


. 097 


2. 0 


. 244061 


. 243057 


. 100 


2. 5 


. 229383 


. 228344 


. 104 


3. 0 


.212563 


.211507 


. 106 


3. 5 


. 194099 


. 193072 


. 103 


4. 0 


. 174496 


. 173578 


. 092 


4, 5 


. 154269 


. 153581 


. 069 


5. 0 


. 133655 


. 133655 


. 030 


5. 5 


. 114105 


. 114382 


. 028 


6. 0 


. 095295 


. 096325 


. 103 


6. 5 


. 078054 


. 079979 


. 192 


7. 0 


. 062821 


. 065697 


. 288 


7. 5 


. 049886 


. 053638 


. 37 5 


8. 0 


. 039362 


. 043751 


. 439 


8. 5 


. 031185 


. 035806 


. 462 


9- 0 


. 025153 


. 029475 


. 432 


9. 5 


. 020963 


. 024413 


. 345 


10. 0 


0. 018243 


0. 020317 


0. 207 



RELATIVE HEATING IN X-Z PLANE- 
9 1 8 • 8M H Z /H Zl=2.0cm 



o 

II 

>- 
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• m m m m • 



o o o o o 



918. 8 mHz FFT RHP 



RELRTIVE HEATING IN X-Z PLANE » 
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o 



o 



o 



918. 8 mHz GQ RHP 
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Table 5. 5 
1500 MHZ 

Comparison of RHPS as Computed by FFT and 



Z = 0. 1 cm Y = 



X 


P FFT 


0 . 0 


0. 769940 


0. 5 


.797484 


1 . 0 


. 782551 


1. 5 


.817937 


2. 0 


. 800998 


2. 5 


. 829700 


3. 0 


. 800512 


3. 5 


. 824841 


4. 0 


. 793202 


4. 5 


. 830136 


5. 0 


. 805427 


5. 5 


. 860066 


6. 0 


. 827499 


6. 5 


. 892359 


7. 0 


. 821465 


7. 5 


. 957599 


8. 0 


. 642953 


8. 5 


. 164369 


9- 0 


. 077163 


9. 5 


. 050834 


10 . 0 


0. 044408 



Z = 1. 0 cm Y = 



0 . 0 


0. 371453 


0. 5 


. 376461 


1 . 0 


. 389622 


1. 5 


. 407001 


2. 0 


. 423560 


2. 5 


. 436171 


3. 0 


. 443940 


3. 5 


. 449051 


4. 0 


. 454231 


4. 5 


. 461729 


5. 0 


. 470444 


5. 5 


. 476507 


6. 0 


. 472567 


6. 5 


. 450246 


7. 0 


. 400007 


7. 5 


. 316297 


8. 0 


0. 212149 



. 0 cm Z i =2.0 cm 

P GAUSS 

0. 791965 
. 795102 
. 803072 
.812381 
. 819285 
. 821610 
. 819920 
. 817286 
. 817557 
. 823257 
> .834203 

. 847315 
. 857974 
. 862831 
. 863583 
. 876068 

.781976 1 

. 155495 
. 083483 
. 061547 
0. 050588 

0. 0 cm Z^ = 2. 0 cm 

0. 375005 
. 379559 
. 391727 
. 407720 
. 423174 
. 435057 
. 442821 
. 448231 
. 453964 
. 461654 
. 470326 
. 475855 
. 471472 
. 448732 
. 398553 
. 315145 
0. 212777 



GQ 



% ERROR 

2 . 202 
0. 238 
2. 052 
0. 556 
1.829 

0. 809 

1. 941 
0. 756 

2. 435 
0 . 688 

2. 878 
1. 275 

3. 048 

2. 953 

4. 212 
8. 153 

3. 902 
0. 887 

0. 632 

1. 07 1 
0. 618 



0. 355 
. 310 
. 210 
. 072 
. 039 
.111 
. 112 
. 082 
. 027 
. 008 
. 012 
. 065 
. 110 
. 151 
. 145 
. 115 
0. 063 
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Table 5. 5 (continued) 



Z = 1. 0 cm Y = 0. 0 cm Z ^ = 2. 0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


8. 5 


0. 128503 


0. 133475 


0. 497 


9- 0 


. 081896 


. 091392 


0. 950 


9- 5 


. 059079 


. 069619 


1. 054 


10. 0 


0. 049324 


0. 056038 


0. 671 




Z = 2. 0 cm Y = 


0. 0 cm Z^ =2.0 cm 




0. 0 


1. 000000 


1. 000000 


0. 000 


0. 5 


1. 023147 


1. 022908 


. 024 


1. 0 


1. 090709 


1. 089751 


. 096 


1. 5 


1. 196503 


1. 194405 


.210 


2. 0 


1. 328885 


1. 325446 


. 344 


2. 5 


1. 470068 


1. 465576 


. 449 


3. 0 


1. 597090 


1. 592548 


. 454 


3. 5 


1. 684959 


1. 682123 


. 284 


4. 0 


1. 711551 


1. 712660 


.111 


4. 5 


1. 662959 


1. 670168 


0. 721 


5. 0 


1. 537580 


1. 552239 


1. 466 


5. 5 


1. 347428 


1. 369424 


2. 200 


6. 0 


1. 115913 


1. 143307 


2. 739 


6. 5 


0. 872545 


0. 901637 


2. 909 


7. 0 


. 646177 


. 672073 


2. 590 


7. 5 


. 458807 


. 476519 


1. 771 


8. 0 


. 321036 


. 327093 


0. 606 


8. 5 


. 230057 


. 224164 


0. 589 


9- 0 


. 172774 


. 158535 


1. 424 


9. 5 


. 133854 


. 117790 


1. 606 


10. 0 


0. 102661 


0. 091728 


1. 093 




Z = 3. 0 cm Y = 


0. 0 cm Z^ = 2. 0 cm 




0. 0 


0. 445169 


0. 443357 


0. 181 


0. 5 


. 454701 


. 452869 


. 183 


1. 0 


. 482492 


. 480589 


. 190 


1. 5 


. 525897 


. 523858 


. 204 


2. 0 


. 579906 


. 577710 


. 220 


2. 5 


. 636850 


. 634624 


. 223 


3. 0 


. 686838 


. 684962 


. 188 


3. 5 


.719186 


.718322 


. 086 


4. 0 


. 724650 


. 725647 


. 100 


4. 5 


0. 697849 


0. 701516 


0. 367 
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Table 5. 5 (continued) 



Z = 3.0 cm Y = 0. 0 cm Z ^ = Z. 0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


5. 0 


0. 639040 


0. 645835 


0. 679 


5. 5 


. 554515 


. 564262 


0. 97 5 


6. 0 


. 455316 


. 467056 


1. 174 


6. 5 


. 354539 


. 366590 


1. 205 


7. 0 


. 263993 


. 274248 


1. 025 


7. 5 


. 191319 


. 1977 68 


0. 645 


8. 0 


. 138623 


. 139991 


. 137 


8. 5 


. 103130 


. 099480 


. 365 


9 0 


. 079439 


. 072423 


. 7 02 


9- 5 


. 062145 


. 054623 


. 752 


10. 0 


0. 047676 


0. 042668 


0. 501 




Z = 4. 0 cm Y = 0. 


0 cm =2.0 cm 




0. 0 


0. 198890 


0. 197535 


0. 135 


0. 5 


. 202728 


. 201401 


. 133 


1. 0 


.213910 


.212657 


. 125 


1. 5 


. 231335 


. 230181 


.115 


2. 0 


. 252897 


. 251864 


. 103 


2. 5 


. 275344 


. 274490 


. 085 


3. 0 


2. 94464 


. 293935 


. 053 


3. 5 


. 305738 


. 305781 


. 004 


4. 0 


. 305365 


. 306291 


. 093 


4. 5 


. 291365 


. 293457 


. 209 


5. 0 


. 264338 


. 267719 


. 338 


5. 5 


. 227505 


. 232026 


. 452 


6 . 0 


. 185905 


. 191090 


. 518 


6. 5 


. 144990 


. 150078 


. 509 


7. 0 


. 109152 


. 113250 


. 410 


7. 5 


. 080784 


. 083091 


. 231 


8. 0 


. 060134 


. 060205 


. 007 


8. 5 


. 045859 


. 043822 


. 204 


9. 0 


. 035893 


. 032519 


. 337 


9- 5 


. 028277 


. 024815 


. 346 


10. 0 


0. 021741 


0. 019488 


0. 225 



RELATIVE HEATING IN X-Z PLANE. 




RELATIVE HEATING IN X-Z PLANE, Y=0 
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Table 5. 6 
2450 MHZ 

Comparison of RHPS as Computed by FFT and GQ 



Z = 0. 1 cm Y = 0. 0 cm Z ^ = 2. 0 cm 



X 


P FFT 


P GAUSS 


% ERROR 


0. 0 


0. 366811 


0. 377543 


1. 07 3 


0. 5 


. 378686 


. 378658 


0. 002 


1. 0 


. 370445 


. 380583 


1. 014 


1. 5 


. 381763 


. 380782 


0. 098 


2. 0 


. 368881 


. 378765 


0. 988 


2. 5 


. 378281 


. 376830 


0. 145 


3. 0 


. 366957 


. 377359 


1. 040 


3. 5 


. 382005 


. 380061 


0. 194 


4. 0 


. 370698 


. 382443 


1. 174 


4. 5 


. 384621 


. 382672 


0. 195 


5. 0 


. 367049 


. 381338 


1. 429 


5. 5 


. 383721 


. 380048 


0. 367 


6. 0 


. 363116 


. 378770 


1. 565 


6. 5 


. 386408 


. 375910 


1. 050 


7. 0 


. 351085 


. 372131 


2. 105 


7. 5 


. 407494 


. 376310 


3. 118 


8. 0 


. 270671 


. 335005 


6. 433 


8. 5 


. 056666 


. 052222 


0. 444 


9. 0 


. 015730 


. 015568 


0. 016 


9- 5 


. 004814 


. 006221 


0. 141 


10. 0 


0. 004538 


0. 004280 


0. 026 




Z = 1. 0 cm Y = 


0. 0 cm Z^ =2.0 cm 




0. 0 


0. 468521 


0. 471441 


0. 292 


0. 5 


. 468444 


4. 71251 


. 281 


1. 0 


. 466835 


. 469570 


. 273 


1. 5 


. 461953 


. 464718 


. 276 


2. 0 


. 454289 


. 457362 


. 307 


2. 5 


. 447671 


. 450967 


. 330 


3. 0 


. 445632 


. 449004 


. 337 


3. 5 


. 448780 


. 451928 


. 315 


4. 0 


. 454143 


.457113 


. 297 


4. 5 


. 458535 


. 461291 


. 276 


5. 0 


. 459520 


. 462096 


. 258 


5. 5 


. 454628 


. 456684 


. 206 


6. 0 


. 437820 


. 439293 


. 147 


6. 5 


. 400307 


. 401186 


. 088 


7. 0 


. 333693 


. 334352 


. 066 


7. 5 


. 238406 


. 238807 


. 040 


8. 0 


0. 135304 


0. 135515 


0. 021 
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Table 5. 


6 (continued) 






Z 


= 1.0 cm Y = 


0. 0 cm Z^ =2.0 cm 




X 




P FFT 


P GAUSS 


% ERROR 


8. 5 




0. 061793 


0. 062571 


0. 078 


9- 0 




. 026235 


. 027940 


. 170 


9- 5 




. 012791 


. 014218 


. 143 


10. 0 




0. 009203 


0. 008925 


0. 028 




Z 


= 2.0 cm Y = 


0. 0 cm Z^ =2.0 cm 




0. 0 




1. 000000 


1. 000000 


0. 000 


0. 5 




0. 974299 


0. 975367 


. 107 


1. 0 




. 907962 


. 911726 


. 376 


1. 5 




. 827885 


. 834780 


. 689 


2. 0 




. 764870 


. 774094 


. 922 


2. 5 




. 7 423 58 


.752317 


• 996 


3. 0 




. 770501 


. 779319 


. 882 


3. 5 




. 844159 


. 850125 


. 597 


4. 0 




0. 942362 


0. 944446 


2. 08 


4. 5 




1. 030636 


1. 029071 


. 156 


5. 0 




1. 069843 


1. 066413 


. 343 


5. 5 




1. 031177 


1. 028739 


. 244 


6. 0 




0. 910033 


0. 911293 


. 126 


6. 5 




. 729405 


. 735520 


.611 


7. 0 




. 529529 


. 539328 


0. 980 


7. 5 




3. 49858 


. 360162 


1. 030 


8. 0 




2. 14015 


. 221086 


0.7 07 


8. 5 




. 1251 18 


. 126845 


. 173 


9- 0 




. 072054 


. 069551 


. 250 


9. 5 




. 040644 


. 037550 


. 309 


10. 0 




0. 021492 


0. 020903 


0. 059 




Z 


= 3. 0 cm Y = 


0. 0 cm Zj = 2. 0 cm 




0. 0 




0. 336450 


0. 336431 


0. 002 


0. 5 




. 328117 


. 328459 


. 034 


1. 0 




. 306858 


. 308100 


. 124 


1. 5 




. 281932 


. 284187 


. 225 


2. 0 




. 263677 


. 266633 


. 296 


2. 5 




. 259506 


. 262592 


. 309 


3. 0 




. 271936 


. 274521 


.2 58 


3. 5 




. 298187 


. 299738 


. 155 


4. 0 




. 330361 


. 330614 


. 02 5 


4. 5 




0. 356574 


0. 355712 


0. 086 
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Table 5. 6 (continued) 



X 



Z = 3. 0 cm Y = 0. 0 cm Z ^ = 2.0 cm 
P FFT P GAUSS 



% ERROR 



5. 0 


0. 364210 


0. 362930 


0. 128 


5. 5 


. 346055 


. 344363 


. 069 


6. 0 


. 299617 


. 300374 


. 076 


6. 5 


. 237352 


. 239797 


. 245 


7. 0 


. 172072 


. 175620 


. 355 


7. 5 


. 115545 


. 118996 


. 345 


8. 0 


. 073442 


. 075554 


.211 


8. 5 


. 045381 


. 045594 


. 021 


9- 0 


. 027663 


. 026528 


. 1 14 


9. 5 


. 016311 


. 015152 


. 116 


10. 0 


0. 008900 


0. 008750 


0. 015 




Z = 4. 0 cm Y = 0. 


0 cm Z^ =2.0 cm 




0. 0 


0. 112552 


0. 112567 


0. 002 


0. 5 


. 109977 


.110109 


. 013 


1. 0 


. 103513 


. 103932 


. 042 


1. 5 


. 096233 


. 096963 


. 07 3 


2. 0 


. 091462 


. 092384 


. 092 


2. 5 


. 091393 


. 092313 


. 092 


3. 0 


. 096492 


. 097207 


. 071 


3. 5 


. 105481 


. 105839 


. 036 


4. 0 


. 115526 


.115471 


. 006 


4. 5 


. 122703 


. 122331 


. 037 


5. 0 


. 123146 


. 122714 


. 043 


5. 5 


. 114733 


. 114582 


. 015 


6. 0 


. 098312 


. 098701 


. 039 


6. 5 


. 077386 


. 078339 


. 095 


7. 0 


. 056350 


. 057606 


. 126 


7. 5 


. 038536 


. 039656 


. 112 


8. 0 


. 025298 


. 025886 


. 059 


8. 5 


. 016291 


. 016220 


. 007 


9- 0 


. 010346 


. 009866 


. 048 


9. 5 


. 006315 


. 005900 


. 042 


10. 0 


0. 003554 


01003536 


0. 002 
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5.4 Analysis of Results 

The different ranges of % ERROR and the number of points lying 
within the various ranges are tabulated in Table 5. 7. Examination of 
the six compared points which lie outside the range of acceptability 
reveals that each of the points has three factors in common 

a. all occur near the edge of the aperture, 

b. all occur when Z = 0. 1 cm, i. e. , Z is small, 

c. all occur in an area where the relative heating pattern 
has a steep slope. 

It was mentioned in section 2.6 that points in the area defined 
by (a) and (b) above may need corrections to account for less than infinite 
integration. These corrections have not been made in the FFT. 

Since only 1. 14% of the points compared fail to meet the 5% 
ERROR criterion and since all of these points are in an area of steep 
slope, the FFT is, therefore, an acceptable computational tool. 



TABLE 5. 7 
Ranges of % ERROR 



Greater 


% ERROR 


But Less 


Number 


Percentage 


Than 




Than 


of Points 


in Range 


10% 






2 


0.38 


5 




10% 


4 


0.76 


4 




5 


8 


1.52 


3 




4 


17 


3.24 


2 




3 


34 


6.48 


1 




2 


71 


15.05 


0 




1 


381 


72.57 



Total 



525 



100 . 00 % 



76 



5, 5 Cost Analysis 

In Chapter 3 it was stated that the fast Fourier transform is an ex- 
tremely efficient method of performing discrete Fourier transforms. In 
this light, two further points are noteworthy in a comparison of the FFT and 
the Guassian guadrature. 

First, all of the computed results in this chapter were performed on 
the CDC 6400 using the RUN (S) compiler. On the basis of central processor 
time only, an average time was taken, over several thousand points, for the 
compution of one point with the following results for the two methods, 

Gaussian guadrature = 1. 9306 sec. /point 5. 5. 1 

Fast Fourier transform = 0. 01499 sec. /point. 5. 5. 2 

This yields a time savings factor (TSF) of 

_ 0.01499 1 

1. 93062 128.8 

which is 7 6% of the theoretical TSF derived in Chapter 3. 

The second point which makes the FFT not only acceptable but also 
desirable is that the FFT computes 4096 points for each x-y plane investigated. 
These points are, for a given Z, in the range 

-16cm < x < 16cm - 16cm < y < l6cm. 

By contrast, the Gaussian guadrature computes only 21 points, for a given 
Z, in the range 



y = o 



0 <x <10cm 
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The advantages of the FFT are, therefore: 

a. ) reasonably comparable accuracy, 

b. ) greater time savings, 

c. ) more complete data, 

d* ) no restrictions on aperture distributions. 



CHAPTER 6 



CONC EUSION 

6. 1 The FFT in Current Applications 

Current development of a safe, efficient, direct contact diathermy 
applicator for use in treatment of musculoskeletal disorders assumes, 
for analysis of the electromagnetic fields, a linear superposition of an 
infinite set of plane wave solutions. The general solutions of the scalar 
wave equation thus obtained develop a complete set of Fourier transform 
pairs which completely specify the electromagnetic field components 
throughout the volume of the human tissue model. 

In the past, the Fourier transforms have been computed by the 
Gaussian quadrature method. Economical use of the Gaussian quadrature 
imposed definite restrictions on the solutions thus obtained: 

a. ) All theoretical and empirical solutions were forced to have 

complete symmetry about the x and y axes thus enabling 
integration in only one quadrant of the u - v plane. 

b. ) The economics of computer time dictated that the relative 

heating be computed only along the positive y-axis in any 
desired z-plane. 

c. ) Extensive corrections were required near the edges of 

the aperture. 

These restrictions obviously limited the generality and completeness 
of any relative heating patterns thus obtained. 
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The fast Fourier transform package developed in this research 
completely circumvents these severe restrictions. Integration is per- 
formed over all four quadrants of the u - v plane, thus eliminating the 
need for symmetry. In addition, the efficiency of the fast Fourier trans- 
form allows the computation of 4096 points in each z-plane rather than 
just a few points along an axis of symmetry. Further, even though no 
corrections are included in the fast Fourier transform package, nearly 
ninety-nine percent of all points calculated agree, to a reasonable 
degree, with answers from the Gaussian quadrature method. The FFT, 
on a point for point basis, accomplishes all of this in two orders of 
magnitude less time. Two minor disadvantages are encountered in the 
use of the fast Fourier transform. Some aliasing is encountered in 
z-planes near the aperture and somewhat more actual computer time 
is required for the computation of each z-plane. Nevertheless, these 
disadvantages are overwhelmingly outweighed by the completeness of 
the relative heating patterns obtained by the FFT. 

6. 2 Future Application of the FFT 

Although the FFT has at this point only been adapted to the 
analysis of relative heating in planar models of human tissue, the 

23 

package developed in this research should be extremely useful to Ho 
in analyzing relative heating in cylindrical models. 

Quantitative studies in the adverse effects of electromagnetic 



radiation will result in the establishment of radiation safety standards. 
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These quantitative studies will require an analysis of the electromagnetic 
fields in both planar and cylindrical models of human tissue. The fast 
Fourier transform package developed in this research can be easily 
adapted to this problem since no symmetry is required in the arrays to 
be transformed. 

These and other problems requiring Fourier analysis in two 
dimensions can be implemented by replacing the SCRPEX and JAMXX sub- 
routines by appropriate subroutines which build the two-dimensional 
arrays to be transformed. 



APPENDIX A 



SUBROUTINE NUMGEN 

21 

This subroutine utilizes that portion of Eguchi's program which 

n k 

generates the weighing factors, W , for each node in the transform 
and calculates the size and beginning and ending values for the DO loops 
used to transform each array. The arrays built by NUMGEN are called 
by FAST in Appendix B. 

If the number of data points is consistent for all uses of this 
package, NUMGEN need only be used the first time the package is used. 
In the first run, the arrays built by NUMGEN can be punched out and 
read in as data in subsequent runs. This will result in slightly less 
central processor time and somewhat more peripheral processor time 
which is desirable on the CDC 6400. Data transmission is via labeled 
common. 
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APPENDIX B 



SUBROUTINE FAST 

21 

Subroutine FAST utilizes that portion of Eguchi’s program 
which calls the subroutine FFT from the bit reversed address arrange- 
ment. FFT is the actual workhorse that performs the Fourier trans- 
formation. Data transmission is via labeled common. 
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APPENDIX C 



SUBROUTINE SCRPEX 

Subroutine SCRPEX calculates the Fourier transform of a 
theoretical aperture distribution for a given u and v. In Chapter 2, 
the final set of equations for the relative heating problem contained a 
| which was the Fourier transform of the aperture distribution at 
z = 0 . This subroutine is entered into the package as a separate 
entity to enable a potential user to change aperture distributions without 
rewriting the JAMXX subroutine. Data transmission is via labeled 
common. 



SUSnOUTI Nh SCHPEX 
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APPENDIX D 



SUBROUTINE JAMXX 

JAMXX computes the appropriate value in the u - v plane for the 
array which is to be transformed. JAMXX requires an input from 
SCRPEX which, for a given u and v, corresponds to theFourier trans- 
form of the aperture distribution. This subroutine by McDougall^ 
incorporates all the theory of Chapter 2 except the corrections and 

the aperture distribution transform. Data transmission is accom- 
plished via labeled common. 
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APPENDIX E 



SAMPLE PROGRAM 

The sample program was used to calculate the relative heating 
patterns in Figure E. 1. It utilizes all of the subroutines in the previous 
appendices in the manner outlined in Figure 4. 1. Figure E. 1 is a 
CALCOMP plot of some of the wealth of information in each transformed 
array. This figure shows the relative heating pattern for various y's 
at the fat-muscle interface for 2450 mHz. It should be mentioned that 
the Gaussian quadrature, for the same problem, obtained only the 
relative heating pattern for y = 0, x_>0. Data transmission with the 
subroutines is via labeled common. 
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